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Using numerical calculations, we compare the collective transition probabilities of many spins in 
random magnetic fields, subject to either frequent projective measurements, frequent phase mod- 
ulations, or a mix of modulations and measurements. For three different distribution functions 
(Gaussian, Lorentzian, and exponential) we consider here, the transition probability under frequent 
modulations is suppressed most if the pulse delay is short and the evolution time is larger than 
a critical value. Furthermore, decoherence freezing (with a transition rate equals to zero) occurs 
when there are frequent phase modulations, while the transition rate only decreases when there 
are frequent measurements and a mix of them, as the pulse delay approaches zero. In the large 
pulse-delay region, however, the transition probabilities under frequent modulations are enhanced 
more than those under either frequent measurements or a mix of modulations and measurements. 

PACS numbers: 03.67.Pp, 75.10.Jm, 03.65.Yz 



I. INTRODUCTION 

Quantum coherence is of key importance in studying 
microscopic and mesoscopic quantum systems and in de- 
veloping quantum devices, including quantum registers 
in quantum computing [lj-l| and spintronic devices @. 
Many methods have been developed to extend the coher- 
ence time of a quantum system, particularly those meth- 
ods using the quantum Zeno effect via either frequent 
measurements or frequent modulations 

Utilizing the quantum Zeno effect, the coherence time 
of a quantum system can be extended sig nificantly 
in nuclear and electron spin systems [10l - fl5j . trapped 
ions [16f 
sates hi 



17 1, ultracold atomic Bose-Einstein conden- 
19 1, and other physical systems [20l - l3l| . Two 



frequent (periodic) control methods are often employed: 
either projective measurements or strong modulations. A 
systematic comparison of these two methods was made 
by Facchi et al. [Hj]: By calculating the transition rates, 
they compared the quantum Zeno/anti-Zeno effect in a 
two-level system via three methods: either (i) frequent 
measurements, (ii) frequent modulations, or (iii) a strong 
coupling to an auxiliary state. They assumed that the 
transition probabilities have always an exponential form 
and thus a well-defined decay rate for short times. 

Zhang et al. [33l - l35| . however, found that some systems 
exhibit decoherence freezing, where the transition proba- 
bility becomes saturated after many modulation periods 
if the pulse delay is short. In these cases, the exponential 
form assumed in Ref. [32[ is violated. Therefore, it is of 
interest to directly compare in the same system the tran- 
sition probabilities (see Fig. [1]), instead of the transition 
rates, under either frequent measurements or frequent 
modulations. 



In this paper, we revisit the study of the transitions of 
many spins in random magnetic fields under either fre- 
quent measurements or frequent modulations. The spins 
will be assumed to be initially in their spin-up state. Us- 
ing exact numerical calculations, we systematically com- 
pare the performance of the three control methods in sup- 
pressing/enhancing the transition probabilities for three 
distributions (Gaussian, Lorentzian, and exponential) of 
the random local fields. 

This paper is organized as follows. In Sec. [TT] we for- 
mulate the spin system's free dynamics and controlled 
dynamics under either (i) frequent measurements, (ii) 
frequent modulations, or (iii) a mix of modulations and 
measurements. We present numerical results in Sec. IIIII 
for the three magnetic field distributions and compare 
in detail the performance of the three control methods 
listed above. Conclusions and discussions are presented 
in Sec. \TV\ 



II. EVOLUTION OF THE SYSTEM 

We consider K spin-1/2 particles in random longitudi- 
nal (z) magnetic fields, but spatially uniform horizontal 
(x) fields |36| : 



H = 



E^k 



fc=l 



K 



ffea 



(1) 



where k is the spin index, and Ukz and Okx are the Pauli 
matrices for the fcth spin. Also, ujk and 2g are the Zee- 
man splitting of the kth spin along the longitudinal and 
transverse direction, respectively. We assume that single- 
spin operations and detection are not accessible, but the 
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FIG. 1: (Color online) Comparison of transition probabilities 
p for free evolution (dashed line), under frequent measure- 
ments (solid line, multiplied by 5), under frequent modula- 
tions (dash-dotted line, multiplied by 100. Presented data 
are for even- number modulation pulses), and under a mix 
of measurements and modulations (dotted line, multiplied by 
100). Here the distribution has a Gaussian form [see Eq. (|18|l ] 
with u} m = and T = 1. The pulse delay is r = 0.2. The tiny 
arrow marks the critical time t c where the transition prob- 
abilities are the same under the modulations and under the 
mix of modulations and measurements. Frequent measure- 
ments and the mix method suppress the transition rate while 
frequent modulations freeze the transition probability. 



ensemble ones are available, which is the case in nu- 
clear spin experiments. We also assume that g is much 
smaller than the typical value of ui k , i.e., !) » j with 
,2/^1/2 



function p(u>), the transition probability becomes 

sin 2 fit 



Pol*) 



g 2 p(Q) Jdoj 
2ng 2 p(Q) t. 



Vl 2 



(3) 



From Eq. ([5]), the transition rate 70 to the spin-down 
state is given by 



To 



dpo 
dt 



g 2 ^ sin2fl k t 



K 



E 

fc=i 



Free evolution 



Note that this transition rate 70 becomes constant if t is 
long enough (but still satisfies t <§C |g| -1 )- Moreover, the 
transition probability po (t) is still small for large enough 
t. In the dense-distribution approximation, the above 
equation for 70 becomes 

70 ~ 2ng 2 J du p(ui)5(u)) 

= 27^(0) (4) 
which is consistent with the result of Eq. ([3]). 



Controlled evolution under frequent 
modulations 



The evolution of two-level systems under many phase- 
modulation contr ol p ulses has been widely investi- 
gated [ll|, [H, HI M, laJSl- The control pulses are 
usually assumed to be hard, collective, and instanta- 
neous [HI], changing the phases of all spins by n. We 
denote such a pulse as a Z pulse. The unitary trans- 
formation of spins due to a Z pulse is described by the 
operator: Z = ® fc (| t)(t I ~ I I) (I |)fc< We have here ne- 
glected the constant i — \f—l which does not affect the 
conclusion. The time evolution operator after N modu- 
lations is mi 



By initially setting all spins in the spin-up state, the 
collective transition probability to the spin-down state 
becomes 



1 K g 2 

= 17 E Sin2 



(2) 



with fl\ = (wfc/2) 2 + <7 2 . If the distribution of uj k is dense, 
we can replace the sum over k with an integral over u> 

f g 2 

Po(t) ~ duj p(u>)— sin 2 Qt, 

where p(ui) is the distribution function and is normalized: 
J du p{u) = 1. In this short time region, <C t <C 
|<?| _1 , with r being the spectrum width, the integrand 
sin 2 Qt/Q 2 sharply peaks at Q « 0. For a flat distribution 



U(t = Nt) = (ZU ) 



N 



—U21 —U22 



N 



(5) 



where 



COS fitT 



z — — — sinSifcT 
2fifc 



U12 = U 2 i = i— sinfi fc T. 

with r being the delay between pulses. After a straight- 
forward simplification, we obtain the controlled-evolution 
of the transition probability from the spin- up state to the 
spin-down state at time t = Nt 



k 



PmodM 



1 V- g 2 ■ 2 o sin 2 iVA fc 

— > — 77 sin \lhT ^ 

K ^.n 2 sin 2 A fc 



k=l 



2 2 sin 2 TV A 

^2 _.__2 x > ( 6 ) 



duj p(uj) --t sin fir — --h— 
sin A 
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where X k is determined by sin X k = 1 — 
(uj k /2n k ) 2 sm 2 n k T. 

In the limit of short time delay r between pulses, r — > 
0, the modulated transition probability p' mod in Eq. ([5]) 
becomes 



P 



1 



mod 



Po(r) 



(7) 



with the transition rate 7 mca s = g 2 T. Compared to the 
modulated case Eq. ([7]), where the transition rate is zero, 
Eq. (|14j) gives a nonzero transition rate 7 mea s, unless r 
is exactly zero. In this sense, as long as the number of 
pulses is large, the transition probability under frequent 
measurements p' mcaB would always exceed that under fre- 
quent modulations (see Fig. [T]). 



with PqM ~ g 2 r 2 , the transition probability in the first 
delay r [35|]. In this limiting case, p' mod becomes inde- 
pendent of the total evolution time and decay freezing 
occurs. 



D. Controlled evolution under a mix of 
modulations and measurements 



C. Controlled evolution under frequent 
measurements 

We now assume that the measurements of the spin 
system are projective and periodic. The effect of such 
a measurement on the spin system is described by the 
projection operator, V — ® k (\ t)(t l)fe- By including 
the free evolution of the system during the measurement 
delay t, we obtain the total evolution operator in a period 
for a single spin as 



V = VU 



U U U12 





(8) 



Note that this evolution is nonunitary because of the 
measurement. It is straightforward to find that for N 
periods the evolution operator is 



V 



N 



TjN 









(9) 



The survival probability of an initially spin-up state for 
the fcth spin becomes 



By combining both frequent modulations and frequent 
measurements, we may utilize the advantages of both 
control methods. Here, a mixed cycle with a period of 2r 
involves a modulation followed by a measurement. The 
nonunitary evolution operator for the cycle becomes 



PU ZU Q 



C/fi - U 12 U 2 i U 12 (Uu - U 22 ) 











(15) 



The total survival probability of an initial spin-up state 
at time t = Nt is 



1 



K 



,2„2 



fe=i ^ k 



N/2 



It is easy to obtain the transition probability 



k=l x K 



N/2 



dui p{oS) 1 



2 2 \ N / 2 

u 9 -in 
— sin i It 

fi 4 



(16) 



PKs - \U^\ 2 = 



1 - 



i A' 



sin 2 ritr 



(10) 



By including all spins, the above Eq. (fTOf becomes 



p s (t = Nt) = 



1 K r 

-Y 

K ^ 

k=l 



1 - 



sin 2 QhT 



N 



(11) 



Under the dense-distribution approximation 



Ps(t) « du p(w) 



1 



^)sin 2 ^ 



(12) 



As a result, the total transition probability away from 
the initial spin-up state becomes 



Pi, 



s (t) = l-p a (t). 



(13) 



In the limit r — > 0, the transition probability in 
Eq. (Tl3|) approaches 



PmeasW ~ 7mcas^ 



(14) 



As seen from Eq. (1T6]) . it is difficult to obtain any ana- 
lytical result in this case without specific information on 
the distribution function p(w). 

In the limit case r — > 0, the transition probability be- 
comes 



Pu 



(17) 



where 7 m i x = (l/2)b 2 g 2 r 3 . When t < t c , the mixed 
transition probability p' mix is smaller than the modulated 
transition probability p' mod - The critical time t c is 

t c = (b 2 r)-\ 

This advantage of the mix method, for short times, is 
shown in Fig. [T] Of course, after many pulses, the tran- 
sition rate of the mix method is nonzero, while that of 
the modulation method is zero. Thus, p' mix > p' 1VLO< \ even- 
tually when t > t c . 

As a summary, in Table [J we list the results for the 
short-r limit for the three control methods . 
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TABLE I: Transition probabilities at t — Nt for three control 
methods (modulations, measurements, and mix) in the limit 
t -> and N > 1. 



,x 10 



Modulations 


Measurements 


Mix 


Pmod = (l/2)ff 2 r 2 


Pmeas = 9* ' ^ '' N 


jC* = (l/2)b 2 g 2 r 4 N 



ILL NUMERICAL RESULTS 

Given an arbitrary t, we have to resort to numerical 
calculations in order to compare the transition probabil- 
ities in Eqs. (J6j, (fT3|) . and (fl6|) . except in the limiting 
case t —¥ 0. Among many forms of distribution func- 
tions p(u), we consider three popular choices: Gaussian, 
Lorentzian, and exponential [3, [35l |44| . 

The Gaussian distribution function used here has the 
form: 



p(u>) = Cexp 



\2 



2T 2 



(18) 



where C is the normalization constant, uj m the peak po- 
sition, and r the spectral width. In the numerical calcu- 
lations, the lower and upper cutoff frequencies used here 
are — w c and w c , respectively, with w c = 100T for the 
three distributions. 

The Lorentzian distribution function used here has the 
standard form: 



C 



The exponential distribution function used is 
p(u>) = C exp 



(19) 



(20) 



There are many ways to compare the performance of 
the three control methods. We employ two methods to 
compare the transition probabilities of the spin system: 
(a) Fix the number N of pulses and varying the pulse 
delay r to investigate the dependence of the performance 
of the control method on the pulse delay; (b) Fix the 
total evolution time t = Nt by varying the number of 
pulses (accordingly the pulse delay r) to investigate the 
dependence on the number of pulses. 



A. Comparison of two-pulse results with different 
pulse delays 

As a starting point, let us compare two-pulse effects via 
either modulations, measurements, or the mix of a mod- 
ulation followed by a measurement. It is easy to find that 
the transition probability subject to two modulations is 



Pmod 



duj p{uj) 



fu 2 

ft 4 



sin Ut 



(21) 




FIG. 2: (Color online) Transition probabilities of a many-spin 
system under two measurements (solid lines) or two mod- 
ulations (dash-dotted lines) for three distribution functions: 
Left column [(a) and (d)], Gaussian; Middle column [(b) and 
(e)], Lorentzian; Right column [(c) and (f)], exponential. The 
peaks of the distributions are chosen as u m = for the top 
row of panels [(a), (b), and (c)] and ui m = 2T for the bot- 
tom row [(d), (e), and (f)], respectively. Dashed lines are the 
transition probabilities for the free evolution. The parameters 
r = 1 and g = 0.001 are used for all cases. Hereafter, time is 
measured in units of T . The transition probabilities under 
two modulations are smaller than those under measurements 
for all rs, if u> m ~ (top panels), and for small ts, if w m = 2F 
(bottom panels). 



and that under two measurements 

-,2 



Pi, 



du> sin 2 ^ T ( 2 



-sm nr 



(22) 



The mix of one modulation and one measurement is ex- 
actly the same as two modulations, p' mix — P^od- 

We plot the two-pulse transition probabilities in Fig. [2] 
for Gaussian, Lorentzian, and exponential distributions. 
For the cases of u) m — 0, the top row of Fig. [2] shows that 
both modulations and measurements suppress the tran- 
sition probability (quantum Zeno effect), compared to 
the free-evolution case. In addition, the transition prob- 
abilities under two modulations p' mod are always smaller 
than those under two measurements p' meas , p' mot i < Pmcas- 
While for the cases of u> m = 2T, both methods (i.e., 
two measurements or two modulations) also suppress the 



transition probability and p' 



< Pn 



if t is small, 



but the two methods enhance the transition probabil- 
ity (quantum anti-Zeno effect) and p' mod > p' mcas , if t is 
large. The cross p' mod — p' mcas occurs around r « 1 in 
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FIG. 3: (Color online) Quantum Zeno effect and anti-Zero 
effect in controlling the transition probabilities with frequent 
measurements (green solid lines), frequent modulations (blue 
dash-dotted lines), and the mix of modulations and measure- 
ments (red dotted lines). The transition probabilities subject 
to control pulses are normalized by dividing the free-evolution 
transition probability for the same time. The horizontal black 
dashed lines in the panels (d, e, f) of the bottom row mark 
the boundary between quantum Zeno and quantum anti-Zeno 
effect. The distribution functions are Gaussian [left column, 
(a) and (d)], Lorentzian [middle column, (b) and (e)], and 
exponential [right column, (c) and (f)]. In the top and the 
bottom row, u> m = and 2F, respectively. The total evolu- 
tion time t = Nt = 10 is fixed. Other parameters are the 
same as in Fig. [2] In the top row, the modulation method 
outperforms the measurement one in suppressing the transi- 
tion probability of the spin system. While in the bottom row, 
the modulation method is worse than the measurement one if 
the number of pulses N is small, but better if N is large. 



B. Comparison of multi-pulse results at fixed 
evolution time 

The transition probabilities of the spin system at a 
fixed time t = Nt = 10 for Gaussian, Lorentzian, and 
exponential distributions are shown in Fig. [3J We nor- 
malize the transition probabilities under measurements, 
modulations, or the mix of modulations and measure- 
ments, by dividing the transition probability of the free 
evolution at the same time. In the top row of Fig. [3J 
where u m = 0, the transition probabilities under fre- 
quent modulations are smaller than those under frequent 
measurements and the probabilities for the mix method 
lie in between, i.e., p' mod < p' miic < p' mcas . In addition, 
Pmod, meas, mix < Po for al l N , showing that frequent mod- 
ulations, frequent measurements, and the mix method all 
suppress the transition probability and only the quantum 
Zeno effect occurs. 

The bottom row of Fig. [3J where u m = 2T, shows a 
more complex phenomenon: 

1. For small N, the quantum anti-Zeno effect appears 
P'mod, meas, mix > Po (i- e -, the enhancement of the 
transition probabilities) for all control methods. 

2. For large N, the quantum Zeno effect 
P'mod, meas, mix < Pa ( Le -> the suppression of 
the transition probabilities) appears for all control 
methods. 

3. For the same value of N, by comparing the perfor- 
mance of the modulation, the measurement, and 
the mix methods, we find p mod > p' m i X > Pmeas 

for 

small N but p' nwd < p' mlx < p' meas for large N. 

4. All three methods intersect around N = 10 (or r = 
1), where p; nod « p' meas « p' mix . 

For extremely small r (large N), which is outside the 
region shown in Fig. [3J the mix method performs better 
than the modulation method. For a given fixed total time 
i, we find that the critical value of r is 



Fig. [U It is interesting that all three distributions show 
consistent and similar results. 

As shown in Fig. [51 for all the cases considered here, the 
transition probabilities increases nonlinearly, for small 
rs, and linearly for large rs. It is proved in Eq. ([3]) 
that the transition probability increases linearly when 
uj^ 1 <t< for the free evolution of the spin system. 
Thus the nonlinearity shown for short ts indicates that 
the quantum Zeno and anti-Zeno effects may occur only 
in this short and nonlinear region. Beyond this short-r 
region, the controlled evolutions are essentially the same 
as the free evolution. 



1 

Correspondingly, the critical number of pulses is 

N c = t/r c = b 2 t 2 . 

We remark here that the above conclusions on the gen- 
eral properties of the quantum Zeno and anti-Zeno effects 
agree with results obtained with other analytical meth- 
ods in certain approximations [32l l38j |. The difference 
between ours and Ref. [H, HI] is that we do not assume 
any specific decay form of the transition probability while 
an exponential form is assumed in those references. In 
fact, an exponential decay form is not always the case, 
especially with the modulation method (see Fig. [TJ). 
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FIG. 4: (Color online) Transition probability of the spin 
system under frequent measurements at the fixed time t = 
10 for different short pulse delays r (black solid line). The 
distributions are Gaussian (left column), Lorentzian (middle 
column), and exponential (right column). Also, we use co m — 
(top row) and u m — 2T (bottom row). The blue dashed 
lines are obtained from Eq. 1)271) . 



small |3jJ |3J, |45|. This relationship shown in Eq. ([55) 
has been verified in Ref. [35j . 

Under the mix control method, the transition proba- 
bility grows linearly 



Pn 



7n 



(26) 



with a well-defined transition rate 7 m ; x = b 2 g 2 T 3 in the 
limit t — > 0. Note that this transition rate depends on 
the distribution function since b 2 = Jdui ui 2 p(ui). 
In the small r limit, it is easy to obtain 



Pn 



P0 



2vrp(0)' 



(27) 



Po 2np(0) ' 

We numerically check the relationship for frequent mea- 
surements by redrawing the large- N results of Fig. [3] and 
using r as the horizontal axis. Figure 2] shows how the 
above limiting results are approached for different distri- 
butions as r approaches zero. 



C. Short r limit at fixed time 

Noticing that all the curves for different distribution 
functions in Fig. [3] give similar tails at large N (small r), 
we next compare how the controlled system approaches 
its limiting case r — > 0. For the free evolution, the transi- 
tion probability of the spin system can be approximated 
as follows if t <C 7(7 1 

Po ~ lot (23) 

with 70 = 2irg 2 p(0) being the transition rate. Similarly, 
the transition probability under frequent measurements 
[see Eq. (13)] at times t < 7~* as is 

Pmcas ~ 7meast- (24) 

The above equation shows that the transition rate 
7mcas = 9 2 t can be defined. Compared to the case of 
free evolution, this rate 7 mea s is independent of the distri- 
bution function and depends linearly on the pulse delay 

T. 

Quite differently, the transition probability [see 
Eq. ([7)] under frequent modulations is frozen for small 

fmod * \9 2 r 2 - (25) 

Obviously, no decay rate can be defined. It is remark- 
able that this transition probability is independent of 
the total evolution time and the distribution function. 
Decoherence is frozen after several control pulses if r is 



IV. CONCLUSION AND DISCUSSION 

In summary, using numerical calculations, we compare 
the transition probabilities of a many-spin system in lo- 
cal random fields (with Gaussian, Lorentzian, or expo- 
nential distributions) under either frequent modulations, 
frequent projective measurements, or the mix of modula- 
tions and measurements. In the small- t region, all three 
control methods suppress the collective transition prob- 
ability of the system. Among the three control methods, 
the modulation one exhibits the largest suppression of 
the transition probability if the total evolution time is 
larger than the critical time, and the transition freezes 
after many modulation pulses. If the time is smaller than 
the critical time, the mix method is the most effective at 
suppressing the transition probabilities. 

In the large r region, all three control methods also 
suppress the transition probability if oj m = 0, but en- 
hance the transition probability if u) m is large. Overall, 
the modulation method changes more drastically the dy- 
namics of the system [25| : The modulation method out- 
performs the other two methods in either suppressing the 
transition in the small r region or enhancing the transi- 
tion in the large r region, provided that the evolution 
time is larger than the critical time t c . 

The modulation-pulse delay is the same in all of our 
calculations. By adopting varying pulse delay, such 
as concatenation or Uhrig's protocol [26|, |27|, we could 
in principle suppress/enhance much more the transition 
probability. We believe that the transition probability 
would deviate more from an exponential decay under 
these more complicated modulation pulse sequences. 
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